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ABSTRACT 

Two dimensional gravitational spectra are derived from simple harmonic 
analysis of range rate tracking data on planetary orbiters. The eccentricity of the 
orbit is arbitrary and results are shown to vary substantially with the aspect angle 
of the tracking line of sight with the orbit plane. The development for arbitrary 
start and stop times (with respect to periapsis) uses modified eccentricity functions 
evaluated by quadrature. 

Simulations with a point-masses model of Venus using tracking data on the 
Pioneer Venus Orbiter show excellent predictions of the average orbiter spectrum 
over one Venus day. The Venus gravitational signal should be above the tracking 
noise level for arc lengths longer than 40° (in true anomaly) about periapsis and 
for terms as high as 55th degree. 

Analysis has been made of tracking residuals from a short arc fit to Mariner 
Mars 9 data over the Hellas Basin (using a complete 6th degree field). Results 
are most consistent with higher residual gravitational power than predicted from 
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Kaula’s rule for Mars. 
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GRAVITATIONAL SPECTRA FROM THE TRACKING OF 
PLANETARY SPACECRAFT IN ECCENTRIC ORBITS 


INTRODUCTION 

The spectrum of a surface harmonic field is a ccmvenient summary of (roughly) the power avaiF 
able at a discrete seris of wavelengths. In the case of the gravitational field at the surface of a planet, 
the wavelength is also a fairly good indicator of the maximum depth of the equivalent spectrum of 
source anomalies (e.g., long wavelength, deep source, see Khan, 1977). Detailed interpretation of the 
spectrum in terms of the depth of density anomalies for the Earth began with Guier and Newton (1965) 
with notable analyses by Allan (1972), Lambeck (1976), and Kaula (1977). Recently Ferrari (1977) 
has applied Allan’s method to the Moon and Mars. My purpose here is merely to show how easily 
the spectrum, in particular a smoothed version of it, can be derived directly from harmonic analyst 
of range rate tracking data to planetary orbiters. This work is an extension of the ideas in Wagner 
and Colombo (1978). Here I consider the ‘fixed direction’ tracking of spacecraft in eccentric 
orbits under arbitrary viewing angles with respect to the orbit plane and periapsis. 

DEVELOPMENT OF THE TRACKING SPECTRUM 

I define the spectrum of the planet’s gravitational field as the potential degree variances: 

£ 

"8 = E + sL)- 

m=0 

where the external potential for a ‘sphericaf planet is given as 

V = — 2^ SgnjSinmX'l , 

^ £=o m=o V r / 

r, <j>, X' being the distance from the center of mass, latitude and longitude respectively (r^ is the 
mean radius, ju the Gaussian gravity constant) and the P^^,^ are fully normalized associated 
Legendre functions (Heiskanen and Moritz, 1967). For convenience, I consider the planet to be 
without rotation and the surface spherical coordinate system to have a pole along the tracking 
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Figure 1. Coordinate System Orientations 


Line of Sight (LOS). The field harmonics Vjjjj are initially given in this system. But, we shall 
see that the invariance of the spectrum in rotation makes this initial choice arbitrary. In what 
follows, primed harmonics refer to the pole of the orbit plane: (a', b')£jjj'. The periapsis of the 
orbit is at tt and co, f and I are the planet centered argument of periapsis, true anomaly and in- 
clination of the orbit to the equator. 


From Hotine (1969), the acceleration of the spacecraft by the gravitation field in the LOS 
direction is given in terms of fully normalized harmonics as: 


LOS = y 


thi^) 

6=0 m=0 \ r / 


6+ 1 


Pj>,n(sin0')[a£n^cosmX' + b^j^sinraX'] , 


( 1 ) 


where. 


Ngjj, = [(2K-l)(fi2 -m2)(2fi + l)“i]i/^and 


( 2 ) 


7 = -^ 


2 



To first order, the orbit of the spacecraft is an ellipse virtually unchanging in a single pats 
of the planet Therefore, (to first order) we want to analyze a time series of the anomalous line 
of sight accelerations along the path of the fixed ellipse. This involves initially a transformation 
of the accelerations to a coordinate system with a pole normal to the orbit plane. FoUowhig 

Kaula (1966, p. 31, 34) (see also Wagner, 1976), I Hnd: 

00 / \ K+l f?(K-tn*evcn) 

LOS = T ^ — j P jn^»(o)[aj„»cosm'(w' + 0 + bg„»sinm'(w' + 01 , (3) 

£-0 ' / tn'-O 


where 


P£m'(o)a£m' 


£ 


L 


iPi, m,(£-m')/2 (I') + P£.m,(£.m')/2(01 



£-m even 
£-m odd 


except 


P£o(o)a£o = 


£ 


E 

m=o 



£*m even 
£-m odd 


for fi even. 


(4) 


(5) 


( 6 ) 


and 

£-m even 

» • 

■®fim ^ £-m odd 

An efficient formula for evaluating the usual inclination functions found in Allan 

(1967a). Here the inclination function^like the associated legendre function^are fully normalized. 


^£m'^°^^£m' ” ^^.m,(£-m')/2^^ ^ “ ^£,m,(£+m')/2^^ \ 

Lafim J 1 


Finally, the anomalous accelerations can be written as a time series by expanding functions of 
r and f in terms of the mean anomaly (M) of the satellite in its orbit of semimajor axis ‘a’ ant» 
eccentricity ‘e*. For a full revolution of data, using [G£pq(e)J Kaula’s eccentricity functions 
(Kaula, 1966, p. 35), I find 

cosm'(tu'+ f)^ ^ r[cosco'm'(G+ + G_)J cosnm 

sinm'(o)' + f)J \sl J \lsinco'm'(G+ + G.)J cosnm 

(7) 

- Isinw'm'(G+ -G_)] sinnm 
+ [cosw'm'(G+-G.)l sinnm 
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where 




( 8 ) 


and the mean anomaly is a linear function of the time (t) from periapsis' M = (2jr/T)t, T being 
the orbit period. 


Let the periodic time series of accelerations for a full revolution be "epresented by the 
Fourier series: 

oo 

LOS = (LOS) + ^ (CL05’„cosnm + SLO^nSinnm) 
n = l 

Then Equations (7), (6), (5) and (4) substituted in (3) gives the Fourier Coefficients of 
Equation (9) as 

°° / r \ K+ 1 £(£-iti even) C ^ even 

CLOS^ = 7 ^ ( - ) coscoW[G^ + GJ ^ [F. + FJ J I 


(9) 


i=o \ ^ / m'=0 

C 

+ sinw 


'm'[G+ + G_] ^ [F_-FjP®”"| 
m = 0 


m = 0 
C-m even 

fi-m odd 


fi-m odd 


(10) 


and 


°p /r \ 8+1 8(8-m even) ^ •\6-meven 

SLOS^ = 7^ ( — ) X) -sinwW[G^-GJ [F- + Fjl 1 

L-l^em J 8 


8=0 \a 


m '=0 
8 


+ COSOJ 


'm'[G^ - G_] [F_ - FJ /’"‘'"’I 

m = 0 V^OmJ 


m = 0 

8 -IT iven 
®8m.J 6-m odd 


8-m odd 


( 11 ) 


where 

“ Fg m 

^ ~ ^, m ,(6 + m ')/2 

except that if m' = 0, ' one F is used in [F_ + F+] . 
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Let the power spectrum of these accelerations be: 

i [CLOSl + SI05JJ . 

What is the expected power spectrum over all possible LOS positions with respect to the planet 
(!', (a) constant)? I have no answer to this question at present But there is an easy answer to 
a complementary ergodic question and that answer appears to be the same. The complementary 
question is: What u the expected spectrum for a given LOS position over all possible gravitational 
fields [C*, S*] whose potential coefficients are chosen randomly, in an uncorrelated manner with: 

where. 

Si = oJ/(2K + 1), 

the mean square coefficient of the planet’s field by degree? 


Squaring Equations (10) and (11) [written for the random field (a*,b*)], and taking the 
expectations (noting that the expectations of cross products of field coefficients are zero), the 
result is: 


~ /r \2(fi+l) «(i~nilwen)/- J£ 

E(I^„(i05)l = 72X;(-) E -^(Gj +G?) E 

K =0 \ ^ / m ’=0 V ',-0 


e 

+ 4(cos2w'm')G+G. ^ 
m=0 



( 12 ) 


where; 

F; = F+, m' 0 
f; = 0, m' = 0 


Note the dependence of this result on the gravitational spectrum (5^). While the reference 
for the field is along the LOS, the surface spectrum is the same with respect to any pole. In 
particular, the in Equation (12) are the same as the field spectrum conventionally referred 
to the planet’s North Pole. 
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NORTH 



Figure 2. Planet, Orbit and LOS Coordinates 


To utilize Equation ( 12 ) for predicting the tracking spectrum, we need to find T and w' 

(the complements of the aspect angles of the view (L) with respect to the orbit plane and the 
periapsis direction) in terms of conventional planet centered coordinates (Fig. 2). Using Figure 
2, I find: 

I' = cos’* [cos I sin + sinlcosL^sinfXp - L;^)] , (13) 

and 

co' = 90® - cos'* •|^[sinn’ 0 sinL 0 + cosir^cosL^cosiL\ - 7 rjy^)]/sinr (14) 

where tt^, it\ are the planet centered latitude and longitude of periapsis, L^, are the cor- 
responding coordinates of the viewing direction (LOS) and Xj, is the planet centered longitude of 
the ascending equator crossing of the orbit. 

The expected spectrum for given aspect angles of the view can be simplified for average or 

more specific viewing conditions. From (12) the expected spectrum averaged over all aspect 

angles (90“-w') with respect to periapsis is simply 

« /r \ 2 (e+i) e(C-^even) i 

ElPj(id5)|<„, = 5^., Y. (Ol *Gl)Y, (t- 

C=0 / ny=0 m = 0 
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(15) 


Oilculations sliow the periapsis aspect has little affect on the aveiage spectrum when full 
realistic fields are considered. However, the spectrum changes by over 100% with viewing aqiect 
to the orbit plane (90"~r). Maximum power occirs when viewing in the orbit plane, minimum 
when the LOS is normal to the plane. This spectrum for the Cross Track (CT) view is especially 
simple to find since for T » 0, the only non zero7is (from Allan, 1967b, p. 1836); 

so that: 


oo /j v2(£+l) «(£Weven) 

-2- Cl 2- (Gj + G2)PL’(<»NL'. <i« 

£>0 ' ® ' 


a generalization of the result for circular orbits in W&gner and Colombo (1978, p. 6). 


The expected spectrum, averaf:ed again over all viewing aspects (!') can be estimated from the 
invariance (Wagner, 1977): 


S (f; + fh = pL'(o)- 


m”0 


Thus in Equation (15), the averaged spectrum is estimated to carry a constant averaged 

, 8(2«-1)2 

= (28+ r* > = 


i(28 + l) 


SO that: 


oo /r\2(£+l) /Ct2C- 1^2 \£(£>m' even) 

— ) E (GJ +G!)Pi„,(o). (17) 

«:^\a/ V3r2£ + 1)/ m'-o 


a similar generalization of the result in Wagner and Colombo (1978, p. 6). 


In general, the anomalous gravitational power spectrum on eccentric orbits contain effects 
at all frequencies from every gravitational term. Here power from a term of degiee £>1 radiates 
to frequericies n greater than £ as well as n less than £> However, as will be seei4 except for low 
degree terms, most of the radiation of power is still downwards. 


7 


THE POWER SPECTRUM FOR INCOMPLETE ARCS OF DATA 

The track spectrum for data 'om circular orbits tends to be uniform over the span of the arc. 
Tlie anomalous accelerations are everywhere at the same distance from the planet. In an eccentric 
orbit, however, the accelerations are smaller away from periapsis, distorting the spectrum. This 
distortion is accounted for by the development for a full revolution of data with the Kaula eccen- 
tricity functions [Equation (7)]. But, for the partial arcs of tracking data which are actually 
taken, these results cannot be used directly. 

One solution io this problem would be to complete a full revolution with artificial data in 
a reasonable manner. A more straightforward approach is to develop the spectrum for data over 
arcs less than a full revolution. For this purpose consider an orbit function defined on the inters 
val Mnjin to M^ax 3). 

The objective functions and their expansions (see Equation (7)] are of the form: 

oo 

<M) = (r/a)'’exp(ivO = ^ X'“*''(e)exp(inM'). (18) 

n= - oo 



Figure 3. Development Interval for Orbit Functions 
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Here, u and v are integers (i - and n are the firequenckss (or wave numbos) of the Fourto* 
expansion in the interval -a to a of the modiAed mean anomaly M'. The notatk)n for the Fourio’ 
coefficients follows that for the Hansen coefficients vdiidi develop the orbit fimctions for full 
revolutions <e.g., Allan, 1967b, p. 1843, 1844). 

To evaluate these modify or partial Hansen coefficients, multiply both sides of (18) by 
particular wave, exp (-iiiM'X and integrate with respect to M* from -a to a. Only the particular 
wave tom remains, yklding: 

1 

X*“'^e) = — / (r/a)“expi(vf - nM')dHr. (19) 

2a 

In general, Mqi^x ^ ^ YfonM* is not an odd function about M' = 0. Thus, for an 

unevenly sampled arc about periapsis, the modified Hansm coefficients are complex numbers: 



1 

r . . 


^x;“*^(e) = 

— 

1 (r/a)“cos(vf - nM )dM 

(20) 


2ir 

Jn 



1 




2x j 

/ (r/a)'*sin(vf - i:M')dM', 
*-ir 

(21) 


and the spectrum of the orbit functions contain twice as many terms as in Equation (7). But 
since most of the useful data on high eccentricity m bits will be fairly close to periapsis on both 
sides, it is convenient to study the simpler case of even sampling centered on periapsis. In this 
case, the modified Hansen coefficients are real only [Equation (20)] as they are for full revolu- 
tions. The development in Equation (7) is then the same with the use of modified Kaula functions 
Gjpq(e) where 

G'jpq(e) = X'g«2p'fq P(e) = ” / (a/O®** cos[(£ - 2p)f - (£ - 2p + q)M'] dM' , (22) 

2x 

and 

hr = . 
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to eccentriaty,' wi5i ^soiwl&te' si^erf, 

^W'thfeir Ta^taitc quadratures de^nd on equal 


to ‘av^ik ^ sidiut^k^ofl^pler^ step, it is i^nvenient to 

change the argument in (22) to eccentric anomaly: 

•E. 


■r\ 10 ^ 2 ’*' ! ••?'-■■ 

each step, it is cc 

. i '. ^ I ,:U ;i1,.' . . It . 


Zp)TWf ’ ll ,1 +^OAV-e)I Vi' 

ifiluotJisq o/(i v;inO”*.^ci o: - -i-: • : . - . •.? , »/■•.! •• •• '. 


(ift) 


Tan E/2 ^ (2 - 2p + q)(E - esinE)\ dE, 


M 


max 


J 




lere 


v>:.- 


^max ^max ~ ®*^^max * 

Ifi '->1 .'■{iJil 5 .‘' “ svl J.':'' '.'i ;'»ji, ". i . ■. •■ . ■ *■ ■' 

SIMULATIONS 

\i) if'iuO-; / ■■■<'- . 5: -■ 5 'w ■ r ; ' , ■ . - i. ^ r . 

Numerous calculations have been made with harmonic and point-mass Earth models which 
^v^fy the track spectrum iJtediCted fof the edcelhtric orbiter. In tliesef analyses, the anomalous 
LOS accelerations along a fixed ellipse are integrated numerically over a given arc length T^i^ in 
time. ' • • ■ ' .••..! 


Ima^ne a set of LOS acceleration data given in this arc. The harmonic analysis of this data 
yields the coefficients in the expansion: ■ ; , .r 

LOS = LOSq + ^ (CL05nCOs27mt/T^+ SZ,05„sin 27mt/Ty^). (24) 

,v ' . ; n=l . . -fi, '■ f;'. . : i- : 


Integrating (24), the LOS velocities (or range^rates) are: 


LdS - LdSfy + LOSqI + (T^ /27rn) CLOSj^ sin 2irnt/T^ 

n=i 


- (Ta /2irn) SZ05„ cos 27rnt/TA • 


(25) 


The secular term LOSqI can be eliminated by subtracting from the velocities a trend line connect- 
ing the first and last points (t = 0, Ty^), which also guarantees continuity of the data there. The 


power spectrum of the resulting (LOS) velocity datr. (n>0) is then related to the power spectrum 
of the scMirce accelerations by: 

ti(LdS) = (2irn/TA)2l^(Lto, (26) 

which also holds for expectations. 

The point-masses (Mascon) model investigated consisted of the central mass (ji) of the Earth 
with 266 randomly chosen concentrations in three layers. The mass magnitudes and locations in 
the shallowest layer (20km depth) are shown in Figure 4 (units: 1.8 x ICT^ii). The locations of 
the masses in the 2 deeper layers are also shown (see also, Wagner and Colombo, 1978, p. 16). 

The gravitational spectrum of this model (Fig. S) is somewhat weaker than Kaula*s rule for the 
Earth (Kaula, 1966, p. 98). 

The fidelity of the spectrum predicted for this model is best revealed in analyses of the funda- 
mental LOS anomalous accelerations themselves. Twenty orbit-arcs were generated with random 
inclinations and nodes for the same (Pioneer-Venus) orbit shape (perigee height: 200km, eccen- 
tricity 0.844, period: 24 hours). Periapsis locations and (infinite distance) viewing aspects (f, cu') 
were also chosen randomly (Fig. 6). All arcs were 180** long in true anomaly (52 minutes) centered 
at periapsis. (The planet did not rotate diuing the arc.) 

For an extreme example, the LOS accelerations in arc 1 , with periapsis directly over a fairly 
strong near-surface concentration, is shown in Figure 7. The power spectrum for this arc (Fig. 8) 
is deficient at low frequencies. The expected spectrum has been calculated for average viewing aspects 
from Equation (17). It gives good predictions for higher frequencies, anticipating the presence of 
some ‘spikey’ behaviour in a typical arc over this model. 

The harmonic analyses of accelerations were done without matching the end points, to con- 
form with the development of the acceler.^tions over a partial revolution. As Figure 9 shows, the 
average power spectrum for the 20 arcs is remarkably well predicted by the theory, even to the 
highest frequencies dominated by the effect of the discontinuity at the end points. 
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Figure 4. Point Mass Magnitudes at 20km Depth for Model Planet* 
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Figure 6. Periapsis Locations and View Aspects for 20 


Acceleration-Data Arcs 
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Figure 7. Anomalous Line of Sight Acceleration — Arc 1 
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MASCON PLANET MODEL 
ORBIT; hp = 200km, e = 0.844, T = 24 HRS. 
ARC LENGTH: 52 MINUTES (CENTERED ON 
PERIAPSIS), SEE FIGS. 4-7. 

• MEASURED FROM DATA IN ARC 1 


FREQUENCY, n: CYCLES/ARC 


Figure 8. Power Spectrum of Simulated Anomalous LOS Accelerations — Arc 1 
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Figure 9. Average Power Spectrum for 20 LOS Acceleration-Data Arcs 
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The analysis of velocity or LOS range rate perturbations requires the data at the ends of the 
arc to be matched. When this is done, predictions of the spectrum are excellent. For example, 

40 similar (180® periapsis centered) arcs of LOS range rate perturbations were analyzed. Here, 
periapsis locations wererandora, but all views were in that direction. The planet model in this case, 
however, was a simple harmonic one consisting of just the complete (0,0) through (2,2)Jerms of 
the Mascon field. Orbit and spectrum computations for this truncated model Were considerably 
faster than with the full field. Also, it was felt the average arc results would converge faster on 
the expected spectrum since severe distortions at periapsis would be avoided (as in arc a^). 

Actually, it appears to take more arcs to determine a good average spectrum for Just these 
low harmonics than for the whole field. Limited local sampling is a notoriously poor way to 
determine global low order harmonics or their average effects. Nevertheless, the results of these 
analyses, with all views in the direction of periapsis, afforded a good test of the predictions for 
various aspect angles. The spectrum from low order terms shows considerably more power with 
along track viewing at periapsis than at any other aspect. The results (in Fig. 10) confirm the 
validity of the aspect dependent parts of the theory [calculated from Equations (12) and (26)]. 

It also shows the considerable amount of power radiated (in an eccentric orbit) to frequencies (n) 
higher than the maximum field degree (fi). 

In the circular orbit case there is no such power radiated upwards except that due to the end 
discontinuity. This latter effect can be minimized by matching higher derivitives (Wagner and 
Colombo, 1978, p. 8, 18). But this adjustment is no longer needed with the present, more gen- 
eral theory for partial revolution data which predicts this high frequency ‘ramp’ distortion. 

SIMULATION OF PIONEER VENUS ORB^ITER MISSION 

In an illustration of results wliich should be forthcoming from tlie Pioneer-Venus mission, 1 
have analyzed the anomalous X<55 velocity residuals in 10 partial revolution arcs over the Mascon 
model planet. The coverage of this data is shown within the dashed lines of Figure 11. Notice 
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II ” 3.96M1 X lO^kiii^/tic^ 

S.m - <0,0) (2^) FROM MA8CON MODEL 

ORBIT: hp-2(IOkm 
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ARCS: 180* LONG (TRUE ANOMALY) 
CENTERED ON FERIAPSI8 

•AVERAGE POWER 8PEC1HUM FOR 40 
RANDOM ARCS VIEWED IN PLANE IN 
PERIAPSIS DIRECTION 


PREDICTION FOR IN PLANE VIEWS 
ALONG-TRACK AT PERIAPSIS 




PREDICTION FOR IN PLANE VIEWS 
\\ IN PERIAPSIS DIRECTION 


PREDICTION FOR NORMAL VIEWS 
OF ORBIT PLANE 


FREQUENCY, n: CY./ARC 


Figure 10. Track rower Spectra for a I-o\v Order 
Gravitational F;ciu Model 
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that only a narrow band north of the equator is covered and the viewing aspect to the orbit 
plane is also limited. But the average spec trum for these arcs is still fairly well predicted (using 
average viewing conditions) as seen in Figure 1 2. 

More detailed estimations for Pioneer Venus have been made with a smooth power spectrum 
of the form: 

$l - AJ«. 

Here 6 = - 2 as in Kaula’s rule for the Earth and A is scaled (from 1 x for the Earth) to the 
size and masses of the terrestrial planets such that interioi stresses are the same as in the Earth 
(Kaula, 1969). Thus, for Mars and the Moon, A should be 6.76 x 10"^ and 36.45 x 1(T®; in fair 
agreement with the actual low degree field (Gapcynski, Tolson, and Michael, 1977; Ananda, 1977). 
The scale rule is: 

A(Terr. Planet) _ /Me\^ (Il\* 

10'^ \Aipj \rej 

so that for Venus it is estimated that A = 1.22 x 10'^. 

Using this smooth spectrum and the elements for Pioneer Venus: Perigee height 200 km, period 
24 hours (eccentricity: 0.844), I have calculated the expected spectrum for two (typical) data 
arcs centered on periapsis (Figs. 13a and 13b). The upper bound line (for average views) includes 
effects of 1st degree harmonics which may be removed with better knowledge of the position of 
the center of mass of Venus. Notice in Figure 1 3a the large variation (up to 200%) of expected 
power with the orbit plane viewing aspect. The variation over periapsis aspects (for the full field) 
is much less significant These figures can be used to estimate the spectrum of the field not in- 
cluded in the orbital data reduction, Le., the truncation power. 

For example, if the 2-way Doppler range rate data has an intrinsic accuracy of 0.03 cm/sec 
for an averaging interval of 1 minute and 26 frequencies are estimated in the 180** (52 minute) 
arc then the ‘white’ noise level for each frequency will be 0.03/(26)*^^ = 0.006 cm/sec/frequency. 
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Figure 12. Average Range Rate Power Spectrum from 
120-Day Pioneer Venus Mission 












POWER (rim); cm par *ac 



Figure 13b. Expected Power Spectrum in LOS Range Rate for a 
Pioneer Venus Orbiter 
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From Figure 13a, this is just at the signal level for 21 cycles. Reading aaoss to lower frequencies 
at this level, this threshold power at 1 cycle is just affected by all the gravitational terms above 
about degree SS. The threshold power at higher frequencies are sensitive to only lower degree 
gravitational effects because the potential at altitude declines strongly with degree. Also, from 
Figures 1 3a and 1 3b, it is noted that appreciable power from low degree harmonics radiate to 
higher frequencies. But for the shorter data-arc at high frequency, this is due more to ieakage’ 
from the end discontinuity than the inherent power of the field. 

THE EFFECT CF ORBIT DETERMINATION ON THE POWER SPECTRUM 

The development of LOS perturbations here is to a fixed ellipse with mean elements which 
have not accommodated themselves to these fluctuations. But particularly when the tracking does 
not cover a full revolution or more, the determined reference orbit will necessarily adjust itself to 
the limited data thereby reducing the overall power of the gravitational fluctuations. In fact, 
tracking coverage is never sufficient; orbits always adjust to accommodate data. 

If the mean orbit were known, the effect of orbit determination is merely the effect on the 
LOS of the variation between the adjusted and the true mean elements of the trajectory. This 
problem has been studied by Gottlieb (1970) in connection with the early direct use of LOS 
accelerations to model lunar Mascons (see also Phillips, et al., 1978). Since the problem arises 
from the^ifference of two ellipses, it can be expected that the greatest distortion (and reduction) 
will be in the lowest frequencies (0, 1 and 2 cycles/arc). On the other hand, if a good low order 
gravity model is used in the trajectory determination, this distorting effect of orbit error will be 
minimized. 

An extreme example of LOS low frequency power reduction in orbit determination is shown 
in Figures 14 and 15. Figure 14 shows the ATS-6 to GEOS 3 LOS range rate signal from the 
Goddard Earth Model 6 (Smith, et al., 1976) less C 2 o> on a 96 minute arc of an elliptic orbit 
before and after adjustment. The average LOS is inclined 60" to the orbit plane and, in spite of 
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Figure 1 4. Simulated ATS-6/GEOS 3 Line of Sight Range Rate Residuals 






the rather long arc, a large change of the orbit plane (mostly) is required in the orbit adjustment, elim- 
inating most of the gravitational variation at 1 and 2 cycles. The end-matched power spectra for 
these two cases are shown in Figure 1 5 together with the predicted spectrum from GEM 6 (less 
€ 20 )* In fact, the power at all frequencies (to the beginning of ‘leakage’) for the unadjusted orbit 
is well predicted by the theory even though both planet and viewing aspects move. Here, Equa- 
tion (15) was used for a full revolution of data. For the adjusted orbit, only the power at 0 (not 
shown), 1 , 2 and 7 cycles is seriously deficient. 

Another, less severe case is shown in Figures 16a,b. Figure 16a shows 40 minutes of LOS 
range rate residuals from ATS-6 to Apollo-Soyuz (rev. 115 over Africa) due to a field complete 
from (10,0) to (45,45). The geopotential was truncated from Rapp’s (1977) model to illustrate 
the excellent predictions of the spectrum for high degree effects. The spectrum (Fig. 1 6b) is for 
the (simulated) residuals on an orbit adjusted to the signal. Evidently the accommodation is 
slight, to (mainly) high frequency fluctuations, because even the low frequency power is well 
predicted. 

The simplest solution to the distortion of the residual tracking spectrum from orbit determin- 
ation is just to ignore the low frequency residual power. As Figures 13a,b show, for the eccentric 
orbiter there should be enough power from low degree terms (C< 3) radiating to higher frequencies 
(n> 2) to maintain sensitivity for their solution. The solution for higher degree terms (C > 3) 
should be easier since they have maximum sensitivity at higher (undistorted) frequencies (n>3). 

A more sophisticated approach would be to solve directly from the tracking data for the 
orbit and point-mass concentrations along or near the track. The LOS accelerations of these on 
the orbit would yield an undistorted power spectrum. Bui, since the smoothed gravitational 
spectrum is all that is sought with limited data, it is probably sufficient for this purpose to ignore 
the lowest frequencies in the interpretation of the track spectrum even for circular orbits. 
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DATA FOR REV. 116 PASS OVER AFRICA 
AFTER ORBIT ADJUSTMENT - DATA DUE 
ONLY TO EFFECTS OF RAPP (1877) GEO- 
POTENTIAL (10,0) -» (45,46) 

•• 

•• 

• • 



TIME FROM BEGINNING OF PASS: SECONDS 

Figure 16a. Simulated ATS-6/Apollo-Soyuz Line of 
Sight Range Rate Residuals 
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Figure 16b. Power Spectrum for Simulated 
ATS-6/Apollo-Soyuz Residuals 
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INVERSION OF TRACK POWER SPECTRA 

The predicted LOS range rate tracking spectrum for an eccentric orbiter can be written as a 
linear transformation of the gravitational spectrum: 

- f; SJ_j65., (27) 

«-0 

where S^ ^ sensitivity coefficients calculable from Equations (12) and (26) for the appropriate 
orbit and viewing conditions. For the ciicular orbit case this system of equations generates two 
simple, uncoupled triangular matrices in 6^ and There is no influence for £ < n in 

this case. Therefore, for i ircular orbits, given a measured track spectrum, the solution for the 
gravitational spectrum can proceed by simple back substitution in (27) starting at a sufficiently 
high degree and frequency. But for the eccentric orbit, the sensitivity matrix is full. Therefore, 
the solution for the gravitational spectrum, using (27), must proceed by inversion of ever larger 
(non symmetric) matrices until the result stabilizes. 


A great simplification can also be accomplished in the eccentric orbit case if only a few 
parameters of a smooth gravitational spectrum is desired. Then Equation (27) can be formulated 
as a set of condition equations for each frequency of the track spectrum in terms of the smooth 
gravitational spectrum parameters. 


For example, suppose we assume the smooth spectrum in the usual form: 

6g = A£®, £>0 
= C, fi = 0 

Equations (27) then become; 

£(max) 

l)2»;n= 1,2,3,.... 
C=2 


(28) 


(29) 


which are non linear in ‘A’ and ‘B.’ 
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But these equations can be linearized about a ‘starting’ solution Aq, Bq, Cq. If the power 
at more than 3 frequencies are given, the equations can then be solved in a ‘least squares’ sense 
to minimize the squared residual IP„ (measured) - Pj(A, B,C)]^. Again, experimentation will be 
required to find whether the results are stable and the residuals smallest as the maximum degree 
considered becomes large. 


The linearized condition equations for A, B and C are: 

dp2 ap2 ap2 

p2 (measured) = p2(Ao,Bo,Cq) + — AA + — AB+ — AC, 


aA, 


aB, 


acn 


(30) 


where: 


ap; , 

" = 2Ao X s2,c(c - 1) 


2Bc 


9A, 


e*2 


gp2 e(max) 

E sj,sce- ip«Loge(e-i) 


3B, 


K=2 


and: 


3P^ 


ac 




(A, B, C) = (Aq + AA, Bq + AB, Cq + AC) 


Of course, if the planets mass is well enough known ‘C’ can be set to zero. I retain it here 
because the Mascon model I use is not adjusted to give 0^0= 0. Two great conveniences of the 
smooth spectrum equations are that: 

(1) we can formulate them to suit the measurements, and 

(2) we can edit or down weight those measurements which have known distortions, for 
example, due to ‘leakage’ at high frequencies, or orbit determination at low frequencies. 


As an example of the numerical solution for a simple gravitational spectrum from a measured 
LOS power spectrum, I have analyzed 10 range rate arcs with respect to the (0,0) through (2,2) 
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model implied by the Mascon planet. The orbits for these arcs were all of high eccentricity, with 
only fair coverage over the planet and viewed from different aspects. Equations (17) and (26) 
were used to calculate the sensitivity coefficients and expected spectrum (Table 1). In the solution, 
only the measured power at n= 1,2, 3 cycles per arc was used. Therefore, the solution for A, B 
and C [iteratively from (30)] was determined without residua! error (starting from 10"^, -2.0 and 
10"^ respectively). The solution took 4 iterations to converge and the result is consistent with the 
discrepancy between the measured and predicted track spectrum; the measured being too high at 
low frequency and too low at high frequency. 

INTERPRETATION OF LOS SPECTRUM: MARS- MARINER 9 

The only planetary spacecraft in highly eccentric orbits for which suitably long arcs of track- 
ing exist (prior to Pioneer Venus) are the Mars-Mariner and Viking orbiters (e.g., Gapcynski, et al., 
1977; Sjogren, et al., 1975). The tracking of these orbiters has resolved gravitational features as 
small as 1200 km (9th degree terms). These features show broad correlation with the topography 
of Mars. From this correlation, Ferrari (1977) has found that over large areas (low degree terms) 
the crust of Mars is significantly stronger (supports significantly higher loads) than that of the 
Earth. This finding can also be inferred indirectly from the low degree gravity spectrum of Mars 
which is significantly stronger than Kaula’s scaled rule. (The assumption must be made, of course, 
that these terms arise from the crust.) It is of great interest to know whether the excess gravi- 
tational power holds for high degree terms as well because if it doesn’t, Mars wculd have a rela- 
tively thinner less developed crust than the Earth. On the contrary, the powerful, broad topography 
of Mars argues for a relatively thicker, more developed (younger?) crust. 

Published residual LOS tracking data (Gapcynski, et al., 1977, p. 4326) has been analj'zed 
from a 194 minute Mars-Mariner 9 pass (rev. 354) over the Hellas Basin (Fig. 17). The range rate 
residual data from Gapcynski, et al., appears to show a significant high frequency signal. The 
residual measurements are referred to an orbit fit to the (one minute averaged) 2 way Doppler 
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Figure 17. Range Rate Residuals: Mars Mariner (9) 





signal with a complete 6th degree field. This field, in turn, used the tracking data in this pass so 
that we should expect to see the effects of aliasing in the residual power spectrum. 

To achieve Figure 1 7, I hand-filtered the original residual points. From the original data the 
noise level appeared to be about 0.03 cm/sec. (I used 6.6cm/sec/Hz to convert the S band Doppler 
(2300MHz carrier) to range rate.) This seemed to be reasonable for 1 minute averaged data Har- 
monic analysis was then performed on a ‘signal’ which was linear (a series of lines) between the 
turning points (highs and lows) of the filtered data. The harmonic noise characteristic of this 
filter was easily determined by averaging 10 analyses of the linear filter with random (a = 

0.03 cm/sec) data at the turning points. (Both signal and noise data were ‘end matched’ before 
being analyzed.) 

The residual spectra of ‘signal’ and noise for this arc are shown in Figure 18 along with the 
expected spectrum for this orbit (e = 0.60, height of periapsis = 1660km, average viewing con- 
ditions) from Kaula’s scaled rule (£> 6). There are 3 parts to this figure, each important to the 
interpretation: 

(1) The spectrum of the ‘signal’ (which includes an unknown amount of noise) declines 
from 1 to 6 cycles then rises abruptly to a maximum at 9 cycles and generally declines 
after that, slowly to about 25 cycles then rapidly. 

(2) The expected spectrum has a maximum at 1 cycle declining slowly to 7 cycles, then is 
almost flat to about 15 cycles where it precipitously loses almost all power. The flat 
portion is due to the ‘sudden’ introduction of gravitational power at £> 6. The rapid 
fall off at higher frequencies is due to the high altitude of the satellite. 

(3) The estimated noise spectrum (shown for an overall level of 0.06 cm/sec) is flat to about 
15 cycles falling steadily after that. Note that 42 turning points are involved, for a max- 
imum estimation of aoout 21 frequencies if the data were etiually spaced in time. The 
declining portion of both the noise and ‘signal’ spectra is clearly a characteristic of the 
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Figure 1 8. LOS Range Rate Residual Power Spectrum for 
Rev. 354, Mars Mariner 9 
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filter. The vertical arrows (in Fig. 14) give upper bounds (I arc in 10) for the fluctuation 
of power among the 10 ‘random’ noise data sets. Roughly speaking, 1 frequency out of 
10 at an average random noise level of 0.06 cm/sec will reach the power level at the end 
of the vertical arrows. 

Gearly, the frequencies beyond 25 cycles are the result of noise which calibrates at some 
level above O.f 6 cm/sec. This result is somewhat surprising unless the data is sampled at, not 
averaged over 1 minute intervals (or is otherwise randomly disturbed). On the other side of the 
spectrum, the low end (n < 6 cycles) is just as clearly due to the gravita*^- lal spectrum above 
6th degree of approximately the level as Kaula’s scaled rule. The fact that the power at 1 cycle 
seems relatively low compared to the others may be the result of orbit adjustment to this data. 

The fact that this low frequency part of the signal is everywhere less than expected may be 
merely due to .aliasing in the low degree gravitational harmonics fit to the Doppler tracking in this 
arc. The field solution has undoubtedly absorbed some of the low frequency information properly 
belonging to higher degree terms. Indeed it is more likely that the gravitational spectrum above 
degree 6:(for this arc) is stronger than Kaula’s rule and the actual aliasing below 7 cycles is more 
severe stUl. This would permit the ‘signal’ highs at 9 and 19 cycles to be explained by the field 
at an average noise level only a little higher than 0.06cm/sec. For these features to be explained 
as noise would require a level of 0.09 cm/sec which would virtually wipe out all of the apparently 
significant low frequency signal. In particular, the high at 9 cycles appears most likely to be a 
gravitational feature since it is the persisting fluctuation which reaches it’s maximum effect at 
periapsis. As a confirmation of the above interpretation, I note that Sjogren, et al., (1975, p. 

2899, 2900) report this data is ndeed sampled (not averaged) at 1 per minute with an average 
noise level of 0.07 cm/sec. 

SUMMARY AND CONCLUSIONS 

The line of sight track spectrum (acceleration: and range rate) due to the gravitational field 
has been simply related t<^ the source spectrum for general planetary orbiters under a variety of 
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viewing conditions. WhOe the theory requires the line of ^ht orientation, the field and the orbit 
to be fixed during a pass, simulations under more relaxed conditions still show excellent results 
(see also Wagner and Colombo, 1978). The expected anomalous gravitational spectrum for the 
Pioneer- Venus orbiter has been worked out in detail for two cases of periapsis centered data. 
Harmonic terms as liigh as degree 55 may be visible in this data. 

It appears likely that the power in the very lowest frequencies of the track spectrum will be 
reduced in the accommodation of the orbit to the data. This effect may either be estimated by 
mass modeling (or error analysis) or the low frequency terms may simply be ignored in the inter- 
pretation of the overall spectrum. 

Interpretation of a residual track spectrum for a Mar»-Mariner9 pass ove” the Hellas Basin 
indicates gravitational features there stronger than Kaula’s (equal-Earth-stress) rule for Mars. 
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